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_ ' Abstract. We have constructed new formulae for generation of solutions for the nonlinear 

' , heat equation and for the Burgers equation that are based on linearizing nonlocal trans- 

r T | ' formations and on nonlocal symmetries of linear equations. Found nonlocal symmetries 

, and formulae of nonlocal nonlinear superposition of solutions of these equations were used 
then for construction of chains of exact solutions. Linearization by means of the Legendre 

— , transformations of a second-order PDE with three independent variables allowed to obtain 

^ ' nonlocal superposition formulae for solutions of this equation, and to generate new solutions 

£^ , from group invariant solutions of a linear equation. 

■ Key words: Lie classical symmetry; nonlocal symmetries; formulae for generation of solu- 
| tions; nonlinear superposition principle 
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■ 1 Introduction 

-t— > ■ 

Nonlocal symmetries of nonlinear equations of mathematical physics stay in the center of at- 
tention of many authors [U (2J [3l HI [5]. Methods used for investigation of differential equa- 
J> . tions (DEs) include application of a wide spectrum of nonclassical and nonlocal symmetries of 

des gala ei ei maim. 

5_j ■ Knowing the symmetries we can construct exact solutions for the equations and then proceed 

with their generating, as well as describe sets of conserved quantities, reduce the initial equation 
to the equations with smaller number of variables and solve other problems related to these 
equations. 

The method of studying of the group invariance properties of the DEs, created by S. Lie, 
has been further developed by G. Birkhoff, L.V. Ovsyannikov, P.J. Olver, W.I. Fushchych, 
N.H. Ibragimov, G. Bluman, J.D. Cole, P. Winternitz and many others. 

In the 1970s V.A. Fok [16] found the symmetry of hydrogen atom in Coulomb field, which, 
as became clear later, could not be found by the Lie method. Then other similar facts started to 
arise, and it required an adequate explanation. Resumption of interest to Backlund transforma- 
tions [T7j in the 1950s and discovery of the Miura transformation allowed explaining of existence 
of an infinite set of conservation laws for Korteweg-de Vries equation, and active period of 
finding of further generalizations of S. Lie theory had started. 

One of such generalizations, conditional symmetries of DEs [23] , that are referred the li- 
terature also as nonclassical or weak symmetries [71 [12], may be found by adding of some 
condition to the equation with this conditional symmetry being the symmetry of the resulting 
system. Usually this condition is an equation of a surface invariant under infinitesimal operator 
admitted by the equation under investigation [131 E] • 
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Theory of Lie-Backlund group transformations was developed by R.L. Anderson and N.H. Ib- 
ragimov in the 1970s [SJ [TOl \TT\ [To] . Such transformations forming a one-parametrical group 
depend on derivatives 

Xi = f (x, «,«,... ,u;e), u a = g a {x,u,u, ...,u;e). (1) 

lk lk 

In these formulae and elsewhere x = (xq set of independent variables, and 

u = {u l ,u 2 , . . . ,u m ) is a set of dependent variables. The corresponding infinitesimal operator 
has the following form: 

X = £ l (x, u,u, . . . , u)di + rf{x } u, u, . . . , u)d u <>, i = 0, 1, . . . , n, a = 1,2, . . . , m. 

1 k 1 k 

Lie-Backlund symmetries defined in such a way present the same concept as generalized vector 
fields in P.J. Olver's terminology [TS]. Dependence of expressions in (TTJ) on derivatives is a reason 
for regarding such symmetry as nonlocal. Note that reversion of the transformation ([I]) requires 
integration procedure. 

Further in the paper we use the following standard notations: 

d^u, n = 0,1,2, ... ,n - 1, u = {u^}, 



dx» ~ " ' ^ j 

= d^d u u, fx, u = 0,1,2, ...,n-l, it = 



Another direction in investigation of nonlocal symmetries is based on representation of the 
equation 

Fi(x, u,u, . . . ,u) = 
l k 

in the form of a conservation law 

Fj(x, u, u, . . . , u) = Di4> l (x, u,u, . . . ,u) + D2(p 2 (x, u,u, . . . ,u) = 

1 k 1 r 1 r 

and on introduction of a new potential variable by the relation 
D2V = (^) 1 (x, u,u, . . . , 

1 r 

It allows studying of a "potential symmetry" of the equation as a classical symmetry of sys- 
tem [JJE] 

Fi(x, u, u, . . . , u) = 0, = </ )1 (x, u, u, . . . , u), D\v = 4> 2 (x, u, u, . . . , u), 

1 k 1 r 1 r 

It may be noted that in the described method for investigation of nonlocal symmetries of DE an 
important place belongs to the recursion operator for an equation introduced by P.J. Olver [18J . 

Further development of this approach was obtained in |19[ I2P] via reversion of the recursion 
operator. Nonlocal symmetries of the pseudopotential type in sense of prolongation structures 
of Estabrook and Wahlquist were considered in [21]. It was discovered in the paper [22] that 
the Jacobi identity for characteristics of nonlocal vector fields "appears to fail for the usual 
characteristic computations" that led to the notion called "ghost symmetries" . That also showed 
impossibility of discussion of Lie algebraic properties for general nonlocal symmetries of DEs. 

Utilization of the method of external differential forms by E. Cartan allowed K. Harrison and 
C.J. Papachristou [1] to calculate classical point (internal or geometrical) symmetry by means 
of so-called isovectors generated by a vector field 

V = C(x,u)di + rf(x,u)d u a. 
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So the internal symmetry of DE was determined by the requirement of invariance of an ideal of 
external differential forms of the system 7j that was generated by the equation under action of 
the Lie derivative £y with respect to vector field V 

£vli = bi k j k + K k Jk + JkMi k . 

Nonlocal symmetry of the given equation was studied in [5] with dependence of factors V on 
derivatives included into this scheme. 

Geometrical theory of nonlocal symmetries of DEs has been developed by A.M. Vinogradov 
and I.S. Krasil'shchik [9] and formulated in the language of jet fiber bundles for the functions 
being solutions of the given equation. 

In this paper we continue investigation of nonlocal symmetries of PDEs [341 [24"] by means of 
nonlocal transformations of variables 

B p (x,u,u,...,u;y,v,v,...,v) = 0, 

1 k 1 t 

p = 1, 2, . . . , n + m, u = {u a }, a = l,2,...,m 
or, in a simpler form 

(xi = h i (y,v,i v), 

T • < Ik 

\u a = H a (y,v,v,...,v). 

v lk 
Under the transformation r the equation 

Fi(x, u,u, . . . ,u) = 

1 n 

is transformed into the new equation 
$(y,v,v,...,v) = 

1 m 

of the order m = n + k. Let us present the obtained equation in the form of equality 
Fi(s, «,«,...,«) = XF 2 (y,v,v,...,v) = 0. 

In Is 

Here A is the differential operator of the order r: r + s = n + k (so it is possible to carry out 
factorization of the equation <3?(y, v , v, . . . , v ) = by means of F 2 (y, v , v, . . . , v ) = 0). In this 

1 m Is 

case we say that the equations 

F\(x, u, u, . . . , u) = and F 2 (y,v,v, . . . ,v) = 

In Is 

are connected by the nonlocal transformation r. 

If invariance algebras of these equations have different dimensions, we can put "extra" sym- 
metries of one equation into correspondence with nonlocal symmetry of another equation. We 
use nonlocal transformations of variables r for construction of formulae for generation of solu- 
tions both in the case of nonlocal invariance of the equation 

Fi(x,u,u,.. .,u) = Fi(y,v,v,.. . ,v) = 

1 n In 

and for other equations. If the equation 
F 2 {y,v,v,...,v) = 

1 s 
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is linear and homogeneous, the transformation r allows construction of nonlinear nonlocal su- 
perposition formulae for solutions of the equation 

F\(x, u,u, . . . ,u) = 0. 

1 n 

The approach to investigation of nonlocal symmetries of DEs described above allows treating 
them as a basis for an algorithm enabling construction of new solutions from one or more given 
solutions. 

In the following section we discuss in more detail application of nonlocal symmetries to 
construction of explicit superposition formulae. 

2 On symmetries and superposition principles 
of nonlinear heat equations 

Nonlinear equations of the class 

u t - d x [h(u) + g(u)u x ] = f(x, u) 

were studied in a number of papers, e.g. [12l [TU [25j H3 [271 [28]. 

Nonclassical (conditional) symmetries of nonlinear heat equations of the general form 

u t = u xx + f(u) 

were completely described in [14]. In [12] nonclassical partial symmetries of equations from the 
family 

ut = {g{u)u x ) x + f(u) 

are constructed. These symmetries are actually conditional symmetries of the considered equa- 
tions [51E1I23]. 

Nonlocal symmetries of some equations from the class 
u t - d x [(p(u,u x )} = 

were found by V.V. Pukhnachev |25[ [26] (in particular, the conditions on the function </>, under 
which the transformation maps the equation into itself). The initial equation can be transformed 
by the Lagrange transformation 

£=/ u(y,t)dy+ (j)(u(0,s),u x (0,s))ds, u(£,t) = [u(x,t)] _1 
Jo Jo 

into an equation from the same class 
oj t = [-oj^uj' 1 , -oj~ 3 uj^. 

Further we will review different approaches to superposition principles for solutions of non- 
linear differential equations. 

One of examples of superposition principles for nonlinear differential equations is given by 
the Bianchi permutability theorem for the sine-Gordon equation, which is adduced e.g. in [32] . 
Permutability of auto-Backlund transformations with different parameters allows constructing of 
an infinite family of soliton-like solutions for the sine-Gordon equation by means of extension of 
two one-soliton solutions which correspond to different values of parameters. The same method 
was applied by H.D. Wahlquist and F.B. Estabrook for the Korteveg-de Vries equation [33] . 
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The idea proposed by S.E. Johnes and W.F. Ames [30] was developed in [29]. Two solutions 

(1) (2) 

u and u of the nonlinear equation 
Fix, u,u,. . . u) = 

1 n 

may be used by means of a procedure 

(1) (2) (1) (2) 

A{u , u) = u * u 

(3) 

to construct a new solution u of the same equation: 

A : ( u , u J — > u . 

(Here * denotes a procedure combining these two solutions into a new one.) 

It was noted in [29] that a superposition principle is a symmetry of the system of two copies 
of the same equation 

rp ( « « &\ n T?l (2) (2) (2 \ n fo\ 

F{x, u,u,...,u)=0, F{x, u,u,...,u)=0. (2) 

In In 

They looked for symmetry operators of system ([2]) of the form 

r = rj{u, u)d m (3) 

whereas the representation symmetry operators of a single equation in the classical approach 
are 

r = 

u 

The corresponding algorithm for generation of solutions is implemented by means of calculating 
of the one-parameter invariance group of system ([2]) 

(3) (1), „,(1) (2) 

u = u = j{u , u;a), 

associated with the operator Here a is the group parameter. The algorithm actually repre- 
sents a realization of the solution superposition principle for the equation under consideration. 
The nonlinear heat equations of the class 

uo + di(Ci(u) + C 2 (u)ui) = (4) 

remained an object of interest for many authors [3 [31 [J3] for a long time. In particular, this 
class of equations @ includes the Burgers equation 

uo + uu\ — u\\ = (5) 

and the nonlinear heat equation 

n - di{u~ 2 ui) = 0, (6) 

which admit linearization by nonlocal transformation of variables. It is well known that the 
equation ([5]) is connected with linear heat equation 

vq - v u = (7) 
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by the Cole-Hopf substitution 
u = — 2(lnu)i. 

The nonlocal transformation which provides linearization of equation ([B]) was first found in |39| 
and then re-discovered in [30]. It can be presented in the form [8] 

u(x ,xi) = — -j— r, xi=v(y ,yi), x = y . 

Here x = (xq,xi) and y = (yo,yi) are tuples of old and new independent variables correspond- 
ingly. This fact was utilized for constructing of the nonlocal generating formulae and formulae 
of nonlinear superposition principle for equation ([6]) in [8] and for ([5]) in this paper. These 
generating formulae are constructed with symmetry operators of the linear heat equation ([7]). 

The invariance algebra of the linear heat equation (|7|) was first calculated by S. Lie and is 
adduced in the standard texts e.g. by Ovsiannikov [37] and Olver [18]. A basis of this algebra 
consists of the operators: 

P = d , Pi = 0i, I = vd v , D = 2y d + yidi, 

G = y Q di-\yivd v , F = y^do - yoVi - \{y{ + 2y )vd v , S b = b(y , yi)d v , (8) 

where b = b(yo,yi) is an arbitrary solution of ([7]), i.e. 6o = &n- 

The operators P\ and G were applied in [8] to construction of generating formulae for equa- 
tion ([5]). The operators Pq and Pi were used for the same purposes for equation ([6]). 

In the present paper new formulae for generating of solutions for equations © and ([6]) were 
found. We also constructed new solutions of equation ([5]) and ([6]) with the nonlocal superposition 
formulae 

u (xq, xi) = —Lo\ m( r + r j, —lo\ m( r ) = u , 

(fc) (fc) , ,(fc) s0 , N 

-2flblnr=ui-5(«) 2 ,ife = l,2. (9) 
for the Burgers equation and [8] 

(3) {1 \ (1) ( !\ (2) (1) (1) (1) (2) (2) (2) 

U [X , X\) = U [X ,T)+U (X ,T), U [Xq, T )dT = U [X , T )d T , 



(1) (2) (k) ( %(k) { % (fc) 

Xi=T + T, T = T\ T U U {X ,T), k = 1,2. (10) 



for equation ([6]) . Here r's denote the functional parameters to be excluded from these formulae. 



3 Formulae for generating of solutions 

We can construct new formulae for generating of solutions of linearizable equations by means 
of combining symmetry operator formulae for generating of solutions of linear equations and 
linearizing transformations. 

Namely, let L be a linear differential equation, F a nonlinear differential equation and r 
a nonlocal transformation which linearizes the equation F to the equation L. If L admits a Lie 
symmetry operator X being a vector field on the space of independent and dependent variables 
then the corresponding differential operator Q acting in the space of unknown functions of L 
(i) 

maps any solution v of L to the other solution 



(2) (1) 

v = Q v . 



(11) 
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(2) 



,(2) 



Then the nonlocal mapping r allows us to obtain a new solution u = r v of F from its 

(l) (l) 

known solution u = r v . Here r is treated in a certain way. Therefore, in the case of 
direct utilization of the transformation r for generating of new solutions of F we have to make 

(1) (1) (2) (2) 

several steps u — > v — > v — > u . At the same time, it is possible to derive a formula for 
generating of new solutions of F without involving solutions of linear equations. This way is 
preferable since it allows to avoid cumbersome calculations and to produce solutions in one step. 

Thus, combining the Cole-Hopf substitution linearizing the Burgers equation ([5]) to the linear 
heat equation © and the action of the Lie symmetry operator Pq = do on solutions of ((7|), we 
obtained the following result. 

Theorem 1. The generating formula for equation ([5]) ; which is associated with the Lie symmetry 
operator Pq of ([7]), has the form 



(2) 
U 



u + 



(1) 

Uq 



i(i) , i,m Y < 



(12) 



Applying the formula (|12p iteratively, we construct chains of solutions of the Burgers equa- 
tion ©: 



1. 



-XQ-Xl 



1 



1 + 



1 



x\ + 2xo + k 
1 



1 + 



xi + 2x + k + 2 



1 + 



x x + 2x + k + 4 
here k is an arbitrary constant; 



3. — 



x\ x\ (6xo — x\ ) xi (60xq — 20xoxf + ^i) 
xq xq (2xo — xf) xq (12xq — 12xoxf + xf) 

x x (840xg - 420xgxf + 42x xf - xf ) 
~* x (I20xg - 180xgxf + 30x xf - xf) ~* ' 



1 — 2 tanh(xo + xx) ■ 



13 + 14 tanh(x + x x ) 121 + 122 tanh(x + xx) 
5 + 4 tanh(x + xx) 41 + 40 tanh(x + xx) 



Consider the Lie symmetry operator Pq + Px = 8q + dx of (JT)) and the corresponding operator 
formulae for generating of solutions of the linear heat equation 



(2) 
V 



(d + dx) { v. 



Theorem 2. The generating formula for equation ([5]) ; which is associated with the Lie symmetry 
operator Pq + Px of 0, has the form 



(2) 
U 



(1) 

u 



(1) , (1) 

Uq + Ux 



(1) 
Ux 



i( u Y + u 



(13) 



Similarly to the previous case, chains of solutions of equation ([5|) can be obtained by for- 
mula (USD: 
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2. 
3. 



-XQ-Xl 



1 



1 + 



1 



x\ + 2xq + k 
2 



1 + 



2xi + 4x + 2k + 3 



1 + 



x\ + 2xq + k + 3 



1 + 



4. — 



Xl 
XQ 



2xi + 4x + 2k + 9 
— 6x0X1 + xf — 2x\ xo + 4xo 
xo(2xo — x\ + 2xiXo) 
24xix[j + 48xg - 4xfxg - 48x 2 x 2 1 - 60xiXg + 4xfx + 20xf x - x\ 
xq(8xq — \.x\x\ — 24xiXq — 12x§ + 4x^x + 12x^xq — xf) 



5. — 1 — 2 tanh(xo + xi) ■ 



23 + 22 tanh(x + x x ) 337 + 338 tanh(x + xi) 
7 + 8tanh(x + xi) ' 113 + 112tanh(x + xi) 



We also combine action of the Lie symmetry operator Pq + Pi on solutions of ([7j) with the 
linearizing transformation of equation ([6]) to the linear heat equation ((7J) . 

Theorem 3. The generating formula for equation ([6]), which is associated with the Lie symmetry 
operator Pq + Pi of ([7]) , has the form 



(x ,xi) = ( ( J) 5 



(u T ) -u (u) -u T (u) 



Xi = — u T { u) + ( u ) 
After applying formula Q14|) to the stationary solution 



(14) 



( m(xo,xi) 



of equation © 

(2), 

u (x ,xi, 



-1 



Cixi + C 2 
we obtain the solution 
-1 



Cixi 



Therefore, this solution is fixed under action of (j!4j) . 

The next chosen operator from the Lie algebra ([5]) is D = 2yodo + yi^i. It acts on solutions 

(2) (1) (1) 

of equation (J7J) as « = 2yo^o + 2/1^1- 

Theorem 4. T/ie generating formula for equation ([5)1, which is associated with the Lie symmetry 
operator D of ([5]), /ias i/ie /orm 



(2) (1) , 
It = It + 



ZXQtiO + « + XlUl 



-X (Ul - o( ) 



.''1 



(1) 



(15) 



2\ w J J 2 "* 

The following chains of solutions of equation ([5]) are obtained with the formula ([151) : 
1 „ 4C 2 x - 4Cxi + 4 



1. C^C + 



C 1 



c + 



C 3 xl - 2C 2 x xi + 4Cx + Cx\ - 2x x 



1 + 



1 



xi + 2x 



1 + 



4x + 2xi + l 



4xq + 4xq + 4xqXi + x\ + xi 



3 El _> El 



x 



x 
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, njL . , s 13x + 5xi + 2 + (14x + 4xi + 4) tanh(x + xi' 

1 — 2 tann(xo + xi) — > — 



5xq + xi + (4xq + 2xi) tanh(x + xi) 



Now we construct new solutions of the equation © by application of the superposition 
principle ©. 



1. For the solutions = — and = 1 H — we have 



2 

E l a;n xi &. 1 2 



(3) Cixie ix ° -C 2 ei 2 -x ( 2 xi + 2x ( 2 ) 
u = ~ — . 

asf q 1 

L xq xi 2. _ 

x (Cie ix o - C 2 e~~~(xQ - x$x{)) 
2. If u = 21, ^ = —1 — 2tanh(x + x\) then 

1 5j^o _i_ x l - 

(3) Cixie 4x o + iC2e~ + ~XQ (2sinh(xo + xi) + cosh(xo + xi)) 

L. 5^0 i ^1 — 

^o(— Cie 4x ° + iC2 cosh(x + xi)e~ + ~XQ ) 



(1) r , ( 2 ) 



3. The solutions V = — , V = _ T 2 _ T1 generate the solution 

Xn ' g ^O 1 ° 



(3) Cixie 4 *o - 2C 2 e X0+Xl xl 

U = ~2 ■ 

xo(Cie~^ - C 2 xl + C 2 4 eX0+a;i ) 
In a similar way, we construct new solutions of equation ([6]) with superposition principle (|10j) . 



Thus, from the stationary solutions « = — — and u = we obtain the solution 



(3) 2C ie 2xo 
u 



-1 - Adxxe 2 ^ + VI + 4CiXie 2a; o 

(1) 1 (2) i 

With other two stationary solutions u = — , u = — 77 we get 



(3) 

It 



f LambertTyf-|(a;i+C)e ;E 0-T L Va;o+^ L , 
-X! + 2eV V 2 J 2 J -C 



Note that the constructed solutions are not stationary. 

All obtained solutions can be extended to parametric families of solutions with Lie symmetry 
transformations or by using other known formulae of generating of solutions. For instance, for 
any real value r and any solution v of the Burgers equation dH) 

(2) „ V X 

v = — 2- 



v + r 



is a solution of the same Burgers equation. 

4 Classical and nonlocal symmetries of equation Slid(w /Xiy ) = 

Consider the equation 

UHU22 ~ u\ 2 - (u qu 22 - ul 2 ) - (uooun - ^oi) = 0- ( 16 ) 
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It is an essentially nonlinear second order scalar PDE in three independent variables, which can 
be linearized to the (l+2)-dimensional d'Alembert equation. Its left hand side is the sum of 
algebraic complements to the diagonal elements of the adjoint matrix to the Hesse matrix (u^ v ) 
with the metric tensor (guv) = diag(l, —1, —1) of the Minkowski space R ' 2 . This equation was 
first investigated in [34]. One of the authors of [M] (W.I. Fushchych) suggested the notation 

Slid(*v) := g uu A^ = 0. 

Here A^ v is the algebraic complement to the element of the Hessian matrix. Hereafter the 
indices [i and v run from to 2, and we assume summation over repeated indices. To the best 
of our knowledge, the Lie symmetries of equation (j 16j) had not been studied. We search for the 
maximal Lie invariance algebra A of equation (|16p by the classical technique [181 [37] . 
In view of the infinitesimal invariance criterion, any operator 

X = ^(x,u)dn + 7](x,u)d u , x = (x ,xi,x 2 ), 

from A satisfies the equation 

X^(SM(u^))\ M = 0, 

where X^ is the standard second-order prolongation of X, M is the manifold determined by 
equation (|16p in the second-order jet space. This is a necessary and sufficient condition for 
infinitesimal invariance of equation (|16p . We confine to the manifold M e.g. with the equality 

u 00 u u - Wqi - u\ 2 + u\ 2 

U22 = ■ 

«11 - 1*00 

and then split the obtained expression with respect to the unconstrained derivatives Uuu, (A 4 , v) ^ 
(2,2), and u^. After solving the derived system of determining equations, we find the general 
form of coefficients of operators from A: 

£° = dx + C 2 xi + C 3 x 2 + C 4 , 
= C 2X + C x xi - C 5 x 2 + C 6 , 
f = C 3 x + C 5 xi + Cix 2 + C 7 , 

T] = CgX + CgXi + CiqX 2 + C\\U + C\ 2 . 

Therefore, the maximal Lie invariance algebra A of equation (|16p is 12-dimensional. A basis 
of A consists of the operators 

P = d , Pi = 0i, P2 = d 2 , P 3 = d u , I = ud u , 

J i = x di + xid , J 02 = x d 2 + x 2 d , J u = -x 2 d\ + x x d 2 , 

D = x d + xidi + x 2 d 2 , Qo = x d u , Qi = xid u , Q 2 = x 2 d u . (17) 

The algebra A can be used for finding of exact solutions of equation (I16p . Exhaustive 
investigation of Lie reductions of equation (fTUI) includes the following steps: 1) construction of 
optimal systems of one- and two-dimensional subalgebras of the algebra A and corresponding 
sets of ansatzes of codimensions one and two; 2) reduction of the initial equations with these 
ansatzes to partial differential equations in two independent variables or ordinary differential 
equations; 3) solving of the reduced equations. 

Here we restrict ourselves to particular reductions with respect to one-dimensional subal- 
gebras. Below we list basis operators of these subalgebras, corresponding ansatzes, reduced 
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equations for the new unknown function ip in two independent variables and some exact solu- 
tions of the initial equation, which are found in the following way. 

1) Pi + Qo- 

u = x xi + (p(x , x 2 ), <pl 2 ~ <Poo<P22 + 1 = 0; 

u = x xi + \x\ + \x\ + C\Xq + C 2 . (18) 

2) Pi + Qi: 

u = \x\ + ip(x ,x 2 ), <fil 2 ~ ^00^22 + ^22 - <A)0 = 0; 
u = \x\t ^{xl-xDixl-xl + Ad) 

T Ci In Ud + x 2 2 -xl + ^{xl-xDixl-xl + Ad^j - \{x\ - x\) + C 2 . 

3) J 02 : 

u = <^(u;,xi), uj = xI~x\, ^{<ft\ ~ <Puu<Pn) + 2u<Pu<Puu - <Pu<Pu + ft = 0; 



-|x? ± \^j{xl -x 2 2 ){xl -xl + ACx) 

± Ci In (2d + x§ - ^ + ^(x2-x2)(x2- a; 2 + 4Ci)') - £(xg - x|) + C7 2 



4) J 2 + Jl2 = 

it = v(uji, 

w l(^12 ~~ ¥>ll¥>22) + 4u! 2 (f 2 (f 22 + 4wi<^2¥>12 + 3(^2 = 



1i = (p(lJ0l,UJ2), U>l=X +Xl, UJ 2 = x\ + x\ - Xq, 



u = Ci ^ + x i ~ x o u= g 2 (x| + xf-xg) + Ci 
x + x\ x + xi 

Another way of application of Lie symmetries to construction of exact solutions is generation 
of new solutions from known ones by symmetry transformations. We can reconstruct the Lie 
symmetry group Q of equation (|16p from the algebra A via solving of a set of Cauchy problems. 
As a result, we obtained that the group Q is generated by the following transformations: 

translations with respect to x and it: x' = x„ + a„, u' = u + 0,3, 



scale transformations with respect to x and u: x'^ = e a4 x^, u' 
addition of a linear in x term to u: x' = x^, u' = u + baXn, 



rotations in the plane (x\,x 2 ): 
x' = xq, x\ = x\ cos ci — x 2 sin ci, x' 2 = x\ sinci + x 2 cosci, u' = u, 

Lorentz rotations in the plane (xq,xi): 
x' = xq cosh c 2 + xi sinh C2, x[ = xo sinh c 2 + x\ cosh c 2 , x 2 = x 2 , v! = u, 

Lorentz rotations in the plane (xq,x 2 ): 
x' = xocoshc3 + X2sinhc3, x[ = x±, x' 2 = xosinhc3 + X2COSI1C3, v! = u, 

where ao, • • • , 05, b^, c\, c 2 and C3 are arbitrary constants. 

Parametrical generation of solutions (|18p and (|19p by means of symmetry transformations 
gives, for example, the following new solutions: 



u = (x 2 a + x )xiy / l - a 2 + |(xoa + x 2 ) 2 + \{x 2 a + x ) 2 + Ci(x 2 a + x ) + C 2 , 

C 2 (sj + x\ -x 2 q) + Ci 
xo cosh C3 + X2 sinh C3 + x\ 
Obtained solutions with parameters can be used in construction of further new solutions, for 
instance, by application of the nonlinear superposition principle. 
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5 Nonlocal linearization and the formula of superposition 
of solutions for equation (1161) 

As established in [34] . the (l+2)-dimensional nonlinear equation Slid(u^) = can be reduced 
to a linear one by the Legendre transformation in the space of variables x = (xq,xx,x 2 ): 

u(x ,x 1 ,x 2 ) = y v + yivi + y 2 v 2 - v, 

X ° = ?;0 ' v = v(yo,y 1 ,y 2 ). (20) 

x\ = vx, 

The Legendre transformation is prolonged to the first and second order derivatives in the space 
of variables x in the following way: 

uo = Vo, ux = 2/1, u 2 = y 2 , 

VXXV 22 ~ v\ 2 V m V 22 - V 2 qVx 2 ^00^22 ~ ^20 
^00 — g , U 0X — g ; ""11 — 2 ' 

VX0V2X - vxxv 2 o v 00 v 2 x - V 20 Vxo v 00 vxx - v^ 2 2 

U02 — g , Uqx — , U 22 — , 

where 5 = det(v AtI/ ) 7^ 0. 

Substituting the obtained expressions for v^ v in (fT6|) . we get the linear d'Alembert equation 

□u = n 00 - uxx - u 22 = 0. (21) 

Since the Legendre transformation is an involution, i.e. it coincides with its inverse transforma- 
tion, application of the Legendre transformation to (|2ip results in the initial equation (|16p . 

In [8] a nonlinear superposition formula was derived for any nonlinear differential equation 
which can be reduced by the Legendre transformation (|20|) to a linear homogeneous equation. 
In the case of equation (I16h this formula has the form: 

u\x)=u\r)+u\e), t£(r) = t§(0), x = r + 9, (22) 

where x = (^0,^1,^2); the subscript ^ of a function u denotes differentiation with respect to 
the /i-th argument of u. r = (tq,tx,t 2 ) and 8 = (0o,8x,6 2 ) are tuples of parameter-functions 
depending on x. 

We apply formula (j22[) for construction of new solutions of (I16p from pairs of its known 
solutions. Let us choose the solutions 

u = -^-2Cx{x\ - (x + ax 2 ) 2 - 2(x + ax 2 )) + 2Cx + C 2 , 

(2) / 

u = -2\ -2Cx((ax + axx + x 2 ) 2 + (axx + ax 2 + x ) 2 ) + C 2 . 



They satisfy the condition det(u^ u ) 7^ 0. Therefore, it is possible to use formula (|22p for finding 
of a new solution of (fT6p . We re-denote the arguments of the first and second solutions as r 
and 9 respectively. Replacing by the expression x^ — 9^, we obtain the solution 



u = -J-2Cx(xx - 9x) 2 + 2Cx(a(x 2 - 9 2 ) + x - 9 ) 2 - (a(x 2 - 8 2 ) + x - 9 ) + 2C 



- 2^2Cx((a9 + a9x + 8 2 ) 2 + {a9 x + a9 2 + 8 ) 2 ) + C 2 . 
The functional parameters 9q, 8x and 9 2 are found from the system 

Rx(-(a(x 2 - 6 2 ) +x - 9 ) - 1) = R 2 {9 a 2 + 9 x a 2 + 29 2 a + 9 x a + 9 ), 
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Ri(xi - 0i ) = R2(0 a 2 + 26 ia 2 + 9 2 a 2 + 8 2 a + a), 

Ri(-(a(x 2 - 8 2 ) + x - 6 ) - 1) = R 2 {6 a 2 + B x a 2 + 29 2 a + x a + 0„), 

where 

i?i = (-2Ci(xi - 0i) 2 + 2Ci(a(x 2 - 2 ) + x - 6 ) 2 - (a(x 2 - 2 ) + x - O ) + 2Ci)~S, 
i? 2 = 2(2Ci((a0 o + a0i +0 2 ) 2 + Hi + a02 + 0o) 2 ))^- 

6 Conclusion 

In this work we continued investigation of nonlocal symmetries of PDEs by means of nonlocal 
transformations of variables. Application of nonlocal symmetries to construction of explicit 
superposition formulae and formulae of generation of solutions was discussed. We obtained new 
formulae for generation of solutions for the Burgers equation and the nonlinear heat equation. 
New superposition principles were constructed for them and then used for obtaining of new 
solutions. The formula of nonlinear superposition of solutions for the multidimensional equation 

Slid(ti A1J/ ) = iinu 22 - u\ 2 - (u 00 u 22 - Uq 2 ) - (u 00 uu - Uqi) = 

was applied to construction of a new solution from two known ones. Such algorithms of gene- 
rating of new solutions represent new nonlocal symmetries of nonlinear equations under inves- 
tigation. 

All obtained solutions can be extended to parametric families by means of Lie symmetry 
transformations or by using other formulae for generation of solutions. 

Acknowledgments 

The authors would like to thank the referees for helpful suggestions and comments. 

References 

[1] Anco S.C., Bluman G., Derivation of conservation laws from nonlocal symmetries of differential equations, 
J. Math. Phys. 37 (1996), 2361-2375. 

[2] Moitsheki R.J., Broadbridge P., Edwards M.P., Systematic construction of hidden nonlocal symmetries for 
the inhomogeneous nonlinear diffusion equation, J. Phys. A: Math. Gen. 37 (2004), 8279-8286. 

[3] Galas F., New nonlocal symmetries with pseudopotentials, J. Phys. A: Math. Gen. 25 (1992), L981-L986. 

[4] Papachristou C.J., Harrison K., Nonlocal symmetries and Backlund transformations for the self-dual Yang- 
Mills system, J. Math. Phys. 29 (1988), 238-243. 

[5] Ibragimov N.H. (Editor), CRC Handbook of Lie group analysis of differential equations. Vol. 1. Symmetries, 
exact solutions and conservation laws, CRC Press, 1994. 

[6] Olver P.J., Rosenau P., The construction of special solutions to partial differential equations, Phys. Lett. A 
114 (1986), 107-112. 

[7] Olver P.J., Rosenau P., Group- invariant solutions of differential equations, SIAM J. Appl. Math. 47 (1987), 
263-278. 

[8] Fushchych W.I., Serov M.I., Tychynin V.A., Amerov T.K., On nonlocal symmetry of nonlinear heat equation, 
Proc. Acad, of Set. Ukraine (1992), no. 11, 27-33. 

[9] Krasil'shchik I.S., Vinogradov A.M., Nonlocal trends in the geometry of differential equations: symmetries, 
conservation laws and Backlund transformations, Acta Appl. Math. 15 (1989), 161-209. 

[10] Anderson R.L., Ibragimov N.H., Lie-Backlund transformations in applications, SIAM, Philadelphia, 1979. 

[11] Ibragimov N.H., Transformation groups applied to mathematical physics, Reidel, Dordrecht, 1985. 



14 



V. Tychynin, O. Petrova and O. Tertyshnyk 



[12] Dzamay A.V., Vorob'ev E.M, Infinitesimal weak symmetries of nonlinear differential equations in two inde- 
pendent variables, J. Phys. A: Math. Gen. 27 (1994), 5541-5549. 

[13] Gandarias M.L., Bruzon M.S., Symmetry analysis and solutions for a family of Cahn-Hillard equations, 
Rep. Math. Phys. 46 (2000), 89-97. 

[14] Clarkson P.A., Mansfield E.L., Symmetry reduction and exact solutions of a class of nonlinear heat equations, 
Phys. D 70 (1993), 250-288, [soTv-int/9306002| 

[15] Ibragimov N.H., Anderson R.L., Lie-Backlund tangent transformations, J. Math. Anal. Appl. 59 (1977), 
145-162. 

[16] Fok V.A., Hydrogen atom and non-Euclidean geometry. Preliminary announcement, Zs. Phys. 98 (1935), 
145-154. 

[17] Lowner C.A., A transformation theory of the partial differential equations of gas dynamics, Nat. Advis. 
Comm. Aeronat. Tech. Notes (1950), no. 2065, 56 pages. 

[18] Olver P.J., Applications of Lie groups to differential equations, Springer- Verlag, New York, 1993. 

[19] Guthrie G.A., Recursion operators and nonlocal symmetries, Proc. Roy. Soc. London A 446 (1994), 107-114. 

[20] Guthrie G.A., Hickman M.S., Nonlocal symmetries of the KdV equation, J. Math. Phys. 34 (1993), 193-205. 

[21] Leo M., Leo R.A., Soliani G., Tempesta P., On the relation between Lie symmetries and prolongation 
structures of nonlinear field equations, Progr. Theoret. Phys. 105 (2001), 77-97. 

[22] Olver P.J., Sanders J., Wang J-P., Ghost symmertries, J. Nonlinear Math. Phys. 9 (2002), suppl. 1, 164-172. 

[23] Fushchich W.I., Serov N.I., The conditional symmetry and reduction of the nonlinear heat conduction 
equation, Proc. Acad, of Sci. Ukraine, Ser. A (1990), no. 7, 24-27. 

[24] Tychynin V.A., Nonlocal symmetry and generating solutions for Harry-Dym type equations, J. Phys. A: 
Math. Gen. 27 (1994), 2787-2797. 

[25] Pukhnachev V.V., Nonlocal symmetries in nonlinear heat equations, in Energy Methods in Continuum 
Mechanics (Oviedo, 1994), Kluwer Acad. Publ., Dordrecht, 1996, 75-99. 

[26] Pukhnachev V.V., Exact solutions of the hydrodynamic equations derived from partially invariant solutions, 
J. Appl. Mech. Tech. Phys. 44 (2003), 317-323. 

[27] Sophocleous C, Transformation properties of variable-coefficient Burgers equation, Chaos Solitons Fractals 
20 (2004), 1047-1057. 

[28] Sophocleous C, Further transformation properties of generalized inhomogeneous nonlinear diffusion equa- 
tions with variable coefficients, Phys. A 345 (2005), 457-471. 

[29] Goard J.M., Broadbridge P., Nonlinear superposition principles obtained by Lie symmetry methods, J. Math. 
Anal. Appl. 214 (1997), 633-657. 

[30] Jones S.E., Ames W.F., Nonlinear superposition, J. Math. Anal. Appl. 17 (1967), 484-487. 

[31] Konopelchenko B.G., On the general structure of nonlinear evolution equations and their Backlund trans- 
formations connected with the matrix non-stationary Schrodinger spectral problem, J. Phys. A: Math. Gen. 
15 (1982), 3425-3437. 

[32] Rogers C, Shadwick W.F., Backlund transformations and their applications, Academic Press, New York, 
1982. 

[33] Wahlquist H.D., Estabrook F.B., Backlund transformation for solutions of the Korteweg-de Vries equation, 
Phys. Rev. Lett. 31 (1973), 1386-1390. 

[34] Fushchych W.L, Tychynin V.A., Generating of solutions for nonlinear equations by the Legendre transfor- 
mation, Proc. Acad. Sci. Ukraine (1992), no. 7, 20-25. 

[35] Fairlie D.B., Mulvey J. A., Integrable generalizations of the two-dimensional Born-Infeld equation, 
J. Phys. A: Math. Gen. 27 (1994), 1317-1324. 

[36] Kumei S., Bluman G.W., When nonlinear differential equations are equivalent to linear differential equa- 
tions, SI AM J. Appl. Math. 42 (1982), 1157-1173. 

[37] Ovsiannikov L.V., Group analysis of differential equations, Academic Press, 1982. 

[38] Bluman G.W., Kumei S., Symmetries and differential equations, Springer, Berlin, 1989. 

[39] Storm M.L., Heat conduction in simple metals, J. Appl. Phys. 22 (1951), 940-951. 

[40] Bluman G., Kumei S., On the remarkable nonlinear diffusion equation (d/dx)[a(u+b)~ 2 (du/ dx)]— du/dt=0, 
J. Math. Phys. 21 (1980), 1019-1023. 



